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O-t. Abstract 

CD . 

We study theoretical aspects of the rotating black hole production and evaporation 
in the extra dimension scenarios with TeV scale gravity, within the mass range in 
^ which the higher dimensional Kerr solution provides good description. We evalu- 

ate the production cross section of black holes taking their angular momenta into 
account. We find that it becomes larger than the Schwarzschild radius squared, 
which is conventionally utilized in literature, and our result nicely agrees with the 
recent numerical study by Yoshino and Nambu within a few % error for higher di- 
mensional case. In the same approximation to obtain the above result, we find that 
the production cross section becomes larger for the black hole with larger angular 
momentum. Second, we derive the generalized Teukolsky equation for spin 0, 1/2 
and 1 brane fields in the higher dimensional Kerr geometry and explicitly show 
that it is separable in any dimensions. For five-dimensional (Randall-Sundrum) 
black hole, we obtain analytic formulae for the greybody factors in low frequency 
expansion and we present the power spectra of the Hawking radiation as well as 
their angular dependence. Phenomenological implications of our result are briefly 
sketched. 



1 Introduction 



The fundamental gravitational scale can be lowered down to TeV scale to 
remedy the hierarchy between the Planck and Higgs mass scales in the large 
extra dimension (ADD) scenario pQ (see also ref. [2] for its stringy realiza- 
tion). 1 In warped compactification (Randall-Sundrum) scenario, both of 
them are scaling together along the location in the warped extra dimension, 
leading again to the TeV fundamental scale at our visible brane |4j. When 
nature realizes such a TeV scale gravity scenario, it is predicted that black 
hole production will dominate over the two body scattering well above the 
fundamental scale, with the geometrical cross section of the order of the 
Schwarzschild radius squared (of the black hole mass equal to the center of 
mass (CM) energy of the scattering) [Hj. Following the observation that black 
holes will mainly decay into the standard model fields on the brane rather 
than into the bulk modes collider signatures of black hole production 
and evaporation are studied comprehensively in ref. [JJ and independently 
in ref. [S]. 2 These two pioneering works are applied in a lot of papers of 
the black hole phenomenology of the ultra-high energy cosmic neutrino sig- 
nature [El E] and of the collider signatures [121 EH Clj . (In ref. [THj, it is 
claimed that the black hole production cross section would be exponentially 
suppressed rather than being geometrical; this is later answered by the semi- 
classical argument |16 4 3 and by the correspondence principle applied to the 
production cross sections of black holes and strings [21]. ) 

The black hole phenomenology opens up the fascinating possibility of the 
experimental investigation of the quantum gravity in the following sense. As is 
emphasized in ref. [JJ , the black hole production hides all the shorter distance 
processes than the Planck length scale behind the event horizon and there 
emerges infrared-ultraviolet duality, i.e., the larger the CM energy becomes, 
the better the semi-classical treatment [25 of the resultant black hole is (since 
the Hawking temperature of it becomes lower). In string theory, where its 
non-perturbative definition is not yet established, this kind of situation (du- 
ality) often appears so that one picture is valid in one limit while the other is 
valid in the opposite limit (see e.g. ref. [23] for review and also refs. [23 E])- 
The region of true interest is the intermediate one at which both pictures 

When the number of extra dimensions is two (and hence their size is around mm), 
rather stringent cosmological constraint M > 100 TeV is imposed [3]. 
2 See also ref. [H] for the study before this observation. 

3 Further claim that the classical black hole formation in the two body scattering is 
proved only with zero impact parameter JJj is answered in refs. |18l 1191 12T)] . 

4 We may observe similar correspondence in the power of the exponential suppression 
of the hard scattering cross section [55], following the argument in ref. |23| . 
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break down and non-perturbative formulation of the quantum gravity (or 
string theory) becomes relevant. Given the status of the theoretical develop- 
ment, experimental signature of quantum gravity at this intermediate region 
would be observed as the discrepancy from the semi-classical behavior in the 
black hole picture valid at high energy limit. Therefore in order to investi- 
gate the quantum gravity effect, it is essential to predict this semi-classical 
behavior as precisely as possible. This is the main motivation of our work. 

After the production phase (the "balding" phase), the black holes are 
well described by the higher dimensional Kerr solution [2Zj if the mass of the 
produced black hole (~ the CM energy of the collision) is large enough to 
neglect the brane tension at the horizon and also small enough to neglect the 
topology and the curvature of the extra dimension(s) 7\. Within the LHC 
energy region, the former condition is satisfied (or marginal) and the latter 
is perfect in the ADD scenario Ej while the former is the same as in the 
ADD scenario and the latter is satisfied in the Randall-Sundrum scenario 
(when the horizon radius is smaller than the curvature length scale which 
is one or two order(s) of magnitude larger than the Planck length scale 5 ) 
[7J[TT]. Throughout this paper, we assume that both the two conditions are 
satisfied. 

The black hole emits most of its quanta (and hence loses most of its 
mass and angular momentum) through the Hawking radiation [21] when 
the above "large-enough" (former) condition is satisfied and hence a few 
hot quanta emitted in the final "Planck" phase, which cannot be treated 
semiclassically, does not consist of the main part of the decay product [Zj. 
(Remember that the smaller the black hole becomes, the hotter the Hawk- 
ing radiation.) In most literature the "spin-down" phase of the black hole 
evolution [7j, in which the black hole shed its angular momentum, is simply 
neglected and the Schwarzschild black hole is used from the start relying on 
the four- dimensional result jH2j that the half life for spin down is a few % of 
the black hole lifetime. 6 To improve this "Schwarzschild approximation", 
it is important to estimate the production cross section of the black holes 
with finite angular momenta. In ref. J33], the production cross section of 
rotating black holes is estimated from the quantum mechanical matrix ele- 
ment between the initial two-plane- waves state and the 'black hole state'. In 
this paper, we take more conservative approach based on the (classical) ge- 
ometrical cross section, 7 in the spirit [THj that a classical description should 

5 The refs. [23 ETJ1 ED] considering mainly the application of the AdS-CFT correspon- 
dence also support this view. 

6 The spin-down phase accounts for about 25% of the mass loss in this four-dimensional 
case 

7 See also ref. for earlier heuristic attempt to estimate the rotating black hole 



2 



be more or less valid for the black hole production in order to avoid the 
Voloshin's exponential suppression mentioned above. 8 

The Hawking radiation is determined for each mode by the greybody 
factor PHE], i-e. the absorption probability (by the black hole) of an incom- 
ing wave of the corresponding mode. 9 Unfortunately, the greybody factors 
have been calculated only for brane- and bulk-scalar modes with the Schwarz- 
schild black hole at present [13]. In current black hole phenomenology, the 
Hawking radiation is either not considered (typically in the cosmic neutrino 
signature) or considered with the greybody factor in the geometrical optics 
limit. 10 To study the evaporation of the higher dimensional black hole 
and to progress the phenomenology further, it is prerequisite to obtain the 
greybody factors of the brane fields (which are the main decay modes of the 
black hole as is mentioned above). 

In this paper, we obtain the brane field equations generalizing the Teukol- 
sky's method in four dimensions EH El EE] We show that they are 
separable into radial and angular parts. For the five- dimensional Kerr black 
hole, we find the analytic formulae of the greybody factors within the low 
frequency expansion. 

In section 2, we present the geometrical production cross section of rotat- 
ing black holes with finite angular momenta in the approximation neglecting 
the balding phase. Our result of the largest impact parameter b max for the 
black hole formation turns out to be in good agreement with the numerical 
result by Yoshino and Nambu [2D] • Within the same approximation, we find 
that the (differential) cross section linearly increases with the angular mo- 
mentum for a given black hole mass (~ CM energy). We also estimate the 
production of the exploding black ring and find that it will possibly form 
when there are many extra dimensions. In section 3, we study the Hawk- 
ing radiation from the rotating black hole. First we derive the brane field 
equations for the spin 0, 1/2, and 1 brane fields from the induced metric 
on the brane in the higher dimensional Kerr black hole back ground and 
show that these equations are separable into radial and angular parts for any 
number of extra dimensions. Next we find the analytic expression for the 
greybody factors for brane fields for the rotating five-dimensional (Randall- 
Sundrum) black hole within the low frequency expansions. We present the 

production. 

8 See refs. [HSIESj for the quantum argument which also claim that Voloshin's suppres- 
sion is not applied. 

9 It is first calculated for spin field [S3, then for spin 1 and 2 fields \M EEH1 HT1 EE?] . 
and finally for spin 1/2 field [^TJ [221 f° r four-dimensional Kerr black hole. 

10 See e.g. refs. ^2^] for consideration of the greybody factor in the geometrical optics 
limit for higher dimensional black hole. 
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power spectra as well as their angular dependence applying these greybody 
factors. In section 4, we present a summary and briefly comment on the 
phenomenological implications of our results. 

2 Production of rotating black holes 

First we briefly review the properties of the rotating (4+n)- dimensional black 
hole |2Zj • Since we assume that the large enough condition (explained in the 
Introduction) is satisfied, the charges of the black hole can be neglected; they 
are at most a few coming from the initial two particles. In general, higher 
dimensional black hole may have \_(n + 3)/2j angular momenta. When the 
black hole is produced in the collision of two particles on the brane, where the 
initial state has only single angular momentum (directed in the brane), it is 
sufficient to consider that the only single angular momentum is non-zero [Zj . 
(This implicitly assumes that the balding phase can be neglected, namely 
that the "junk" emissions are negligible and do not change the |_(n + 3)/2j 
angular momenta during this phase.) In the Boyer-Lindquist coordinate, the 
metric for the black hole with single angular momentum takes the following 
form [2Z| 

9 - f 1 - # - si - 2 * ( r2 + ° 2 + ° 2 si - 2 v 

+ 2a sin 2 $^-^-dtdip - ^ri^-dr 2 - E(r, $)d$ 2 
S(r, v) A(r) 

-r 2 cos 2 $dQ n , (1) 

where 

S(r, $) = r 2 + a 2 cos 2 
A(r) = r 2 + a 2 - fxr~ n+l . 

We can see that the horizon occurs when A(r) = 0, i.e. when r = with 

l/(n+l) 

= (l + a 2 )- 1 ^ +1 V 5 , (2) 

where a* = a/r h and the Schwarzschild radius r s = / u 1 ^ n+1 - ) are introduced 
for later convenience. Note that there is only single horizon when n > 1 
(contrary to the four-dimensional Kerr black hole which has inner and outer 
horizons) and its radius is independent of the angular coordinates. We can 



l + a 2 
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obtain the total mass M and angular momentum J of the black hole from 
the metric (|T|) 



M 



(n + 2)A n+2 
16nG 



J 



n + 2 



2 



Ma, 



(3) 



where A n+2 =2 7r( n+3 ^ 2 /r( ! ^) is the area of unit sphere S n+2 and G is the 
(4 + n)-dimensional Newton constant. Therefore we may consider [i and a 
(or r^ 1 and a*) as the normalized mass and angular momentum parameters, 
respectively. We note that there are no upper bound on a when n > 2 nor 
on a* when n > 1, contrary to the four- dimensional case where both a and 
a* are bounded from above. In this paper, we concentrate on the brane 
field equations and hence only the induced metric on the brane is relevant, 
where the last term in eq. (PJ) vanishes and the angular variables $ and (p are 
redefined to take the values < $ < 7r and < <p < 2n. The explicit form is 
given in eq. (joTIj) . 

2.1 Production cross section 

We estimate the production cross section of rotating black holes within the 
classical picture. Let us consider a collision of two massless particles with 
finite impact parameter b and CM energy -y/s = Mj so that each particle has 
energy Mi/2 in the CM frame. The initial angular momentum before collision 
is J, i = bMi/2 (in the CM frame). Suppose that a black hole forms whenever 
the initial two particles (characterized by Mj and J«) can be wrapped inside 
the event horizon of the black hole with the mass M = Mj and angular 
momentum J = Jj (see Fig. Qfor schematic picture), i.e., when 



where rh(M, J) is defined through eqs. (J2J) and Since the right hand side 
is monotonically decreasing function of b, there is maximum value 6 max which 
saturates the inequality (jlj) 



where r s {M) is defined by r s {M) = /i(M) 1 /( n+1 ) and eq. ©. When b = 6 max , 
the rotation parameter a* takes the maximal value (a*) ma x = (n + 2)/2. 

The formula (JSJ) fits the numerical result of 6 max with full consideration 
of the general relativity by Yoshino and Nambu [20] within the accuracy less 



b < 2r h (M, J) =2r h (Mi, bMi/2), 



(4) 




(5) 
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than 1.5% for n > 2 and 6.5% for n = 1 (although it just gives the Schwarz- 
schild radius 6 max = r s (M) for n = which is 24% larger than the numerical 
result HE!): 



n 





1 


2 


3 


4 


5 


6 


7 


-^Numerical [20] 


0.804 


1.04 


1.16 


1.23 


1.28 


1.32 


1.35 


1.37 


R Analytic 


1.00 


1.11 


1.17 


1.22 


1.26 


1.30 


1.33 


1.36 



where R denotes R = fe max /rs(M). 

Our result is obtained in the approximation that we neglect all the effects 
by the junk emissions in the balding phase and hence that the initial CM 
energy Mj and angular momentum Jj become directly the resultant black 
hole mass M = Mi and angular momentum J = Jj. 11 The coincidence 
of our result with the numerical study 20\ suggests that this approximation 
would be actually viable for higher dimensional black hole formation at least 
unless b is very close to 6 max - 12 

Once we neglect the balding phase, hence the junk emission, the initial 
impact parameter b directly leads to the resultant angular momentum of the 
black hole J = bM/2. Since the impact parameter [b,b + db] contributes to 
the cross section 27rbdb, this relation between b and J tells us the (differential) 
production cross section of the black hole with its mass M and its angular 
momentum in [J, J + dJ] 



where 



with 



13 



Jn 



da(M, J) 



Jr\ 



2 « vr ^r(^) 



%-nJdJ/M 2 




femaxM 



(J < J n 
(J > J n 



M 



(n + 2) 



1+ (*±* 



l/(n+l) 



Mr 



n + 2 
n+1 



(6) 



(7) 



Mr 



'(2nT 
8ttG 



l/(n+2) 



(8) 



11 The authors of ref. [201 have found that the irreducible mass of the black hole is 
substantially reduced when b is close to b max and have suggested that balding phase is not 
negligible when b ~ b max . However, the irreducible mass provides the lower bound on the 
final mass of the black hole; at this stage we cannot conclude how much junk energy and 
angular momentum are radiated at the balding phase. 

12 See refs. ^] 03 for estimation of the energy loss during the balding phase for 
the head-on collision (b=0) case obtained from gravitational radiation emitted during the 
infall of a particle into a four dimensional black hole. 

13 See e.g. ref. JHj for different conventions for Mp. 
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The numerical values for j n are summarized in Table H 
Table 1: Numerical values for j n and k r , 



n 





1 


2 


3 


4 


5 


6 


7 


Jn 

k n 

kn/ jn 


0.0398 
0.0159 
0.399 


0.256 
0.125 
0.489 


0.531 
0.228 
0.429 


0.815 
0.251 
0.308 


1.09 
0.214 
0.195 


1.37 
0.155 
0.114 


1.63 
0.101 
0.0619 


1.88 
0.0603 
0.0320 



It is observed that the differential cross section © linearly increases with 
the angular momentum. We expect that this behavior is correct as the first 
approximation, so that the black holes tend to be produced with larger an- 
gular momenta. At the typical LHC energy M/Mp = 5, the value of J max is 
■/max = 2.9, 4.5, . . . , 10, 12 for n = 1, 2, . . . , 6, 7, respectively. This means that 
the semi-classical treatment of the angular momentum becomes increasingly 
valid for large n. 

Integrating the expression © simply gives 



a(M) = 

= Firr s (Mj 2 . 
The factor F is summarized as 



1 + 



n + 2 



-2/(n+l) 



irr s (Mf 



(9) 



n 





1 


2 


3 


4 


5 


6 


7 


-^Numerical |20- 


0.647 


1.084 


1.341 


1.515 


1.642 


1.741 


1.819 


1.883 


-^Analytic 


1.000 


1.231 


1.368 


1.486 


1.592 


1.690 


1.780 


1.863 



This result implies that, apart from the four-dimensional case, we would 
underestimate the production cross section of black holes if we did not take 
the angular momentum into account and that it becomes more significant for 
higher dimensions. We point out that this effect has been often overlooked 
in the literature. 



2.2 Rotating black ring 

In four dimensions, the topology of the event horizon must be homeomorphic 
to two-sphere and there is a uniqueness theorems for static or stationary black 
holes. On the other hand, a higher- dimensional black hole can have various 
nontrivial topology and the uniqueness property of stationary black 



7 



holes fails in five (and probably in higher) dimensions. The typical example 
in five dimensions has been recently given by Emparan and Reall jH] ■ They 
have explicitly provided a solution of the five- dimensional vacuum Einstein 
equation, which represents the stationary rotating black ring (homeomorphic 
to S 1 x S 2 ). In this case, the centrifugal force prevents the black ring from 
collapsing. When the angular momentum is not large enough, the black ring 
will collapse to the Kerr black hole due to the gravitational attraction and 
some effective tension of the ring source. In fact, this five dimensional black 
ring solution has the minimum possible value of the angular momentum given 
by 

/ M \ 3/2 

Jjnhx = &BR \J^) ' ( 10 ) 

where &br = 0.282. On the other hand, we have the upper bound for the 
angular momentum of the black holes produced by particle collisions: 

/ M \ 3/2 

J^ x = Ji{ Wp ) , (11) 

where ji = 0.256. Since these numerical values are of the same order, we 
cannot conclude the possibility of black ring productions at colliders. 14 

Here we consider the possibility of the higher dimensional black ring which 
is homeomorphic to S 1 x S". Corresponding Newtonian situation will be the 
system of a rotating massive circle. They are always unstable in higher 
dimensions; a circle with slow rotation collapses and one with rapid rotation 
explodes toward infinitely large thin circle. In general relativity, we have no 
idea as to the validity of this picture due to the nonlinearity of the Einstein 
equation. We shall discuss in the following the possibility of the black ring 
formation based on Newtonian picture assuming that the nonlinear effects of 
the gravity unchange the qualitative feature. For simplicity, we just consider 
the gravitational attraction and the centrifugal force of the massive circle 
and neglect the effect of tension. Let £, M and J be the radius, the mass 
and the angular momentum of the massive circle. Then we obtain (3 + n)- 
dimensional effective theory with the Newton constant G/27r£ by integrating 
out along the ^-direction. The Schwarzschild radius of the point mass in 
the effective theory is given by 



r 



16ttG M 
(n + l)A n+1 2tt/ 



l/n 



8GM 



{n + l)L4 n+1 



l/n 



(12) 



14 Even if the black rings are produced, they might be unstable due to the existence of 
J m i n and the black string instability. D.I. is indebted to Roberto Emparan for this point. 



s 



Thus we expect the black ring with S 11 -radius t and >S n -radius r. In flat space 
picture, t > r should hold for black ring. This condition gives 



> / ■ 



8GM 



(n + 1) A.+1 



l/(n+l) 



(13) 



On the other hand, the condition that the centrifugal force dominates against 
the gravitational attraction becomes 



j > 2~ (n+3)/2 G 1/2 £~ (n_1)/2 M 3/2 . 



(14) 



This combined with eq. (|13|) gives the minimum value of the angular mo- 
mentum for exploding black ring: 



7 > 7 ■ 



rsj ^mm ,x "n 



M . (n+2)/(n+l) 



Mf 



(15) 



where 



2re 2 +3n + 7 (n+6)(ra-l) 
2 2(n+l) 7]- 4(n + l) 



n 



2(n+l) 



(16) 



The numerical values for k n are presented in Table ^ This result shows that 
Jmm for exploding black rings is one or two order (s) of magnitude smaller 
than J max for collision limit when n is large. Therefore we expect that the 
exploding black rings are possibly produced at colliders if there are many 
extra dimensions, though they will suffer from the black string instability 
when they become sufficiently large thin rings. In the following of this paper, 
we do not follow the evolution of the exploding black ring nor consider the 
radiations from it since this is still at the heuristic stage; we concentrate on 
the Hawking radiations from the higher dimensional Kerr black hole after 
the balding phase. 



3 Radiations from rotating black hole 

In this section, we study the Hawking radiation [21] from the higher di- 
mensional Kerr black hole |2Zj- The Hawking radiation is thermal but not 
strictly black body due to the frequency dependent greybody factor T, which 
is identical to the absorption probability (by the hole) of the corresponding 
mode I31j. The quantity 1 — V for each mode can be computed from the 
solution (to the wave equation of that mode) having no outgoing flux at the 
horizon as the ratio of the incoming and outgoing flux at infinity. 
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It can be shown that the higher dimensional black hole radiates compa- 
rable amount of energy into one brane mode and into one bulk mode (with 
all the Kaluza-Klein tower summed up) ;6j. Typically, the number of degrees 
of freedom is much larger for brane mode than for bulk mode, i.e., tens of 
the standard-model degrees of freedom are living on the brane while there 
are only few degrees of freedom of the graviton (and possibly other fields) 
in the bulk. Therefore the higher dimensional black hole radiates mainly on 
the brane [Hj. For this reason, we concentrate on the greybody factors for 
the brane mode in this paper. 15 



3.1 Brane field equations 

We derive the wave equations of the brane modes using the induced four 
dimensional metric of the (4 + n)-dimensional rotating black hole [27] • The 
wave equations can be understood as generalization of the Teukolsky equa- 
tion [HHl EU El EE] to the higher dimensional Kerr geometry. The derivation 
is shown in Appendix. 

We present the brane field equations for massless spin s field which are 
obtained from the metric (JIJ) with the standard decomposition 

$ s = R s (r)S($)e- wt+m v, (17) 

utilizing the Newman- Penrose formalism [21] 

1 d ( . JS\ 
sin 17- 



sin $ d$ \ d$ 



+ 



[s — atu cos -#) 2 — (scot-# + mcscd) 2 — s(s — 1) + A 



S = 0, 



dr \ dr I 




+ 2maui — a 2 oj 2 — A 



R = 0. 



where 



,[2r+ (n- l)ur- n ]K 

i- — n(n — l)fir~ 



K 



(r 2 + a 2 )u) — ma. 



(20) 



The solution of eq. (fTHj) is called spin- weighted spheroidal harmonics s Si m 
(see e.g. ref. [SIE2]) which reduces to the spin- weighted spherical harmonics 



15 We note that the bulk graviton emission may not be negligible for highly rotating 
black holes since the superradiant emission is more effective for higher spin fields 
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s Yi m (see e.g. ref. jo2]) in the limit aoo <C 1, 

s 5 ;m (aw;i5,(/j) = a YJ m (#, </?) + 0{au) 

where 16 

~(l + m)\(l-m)\2l + l 



(21) 



' -^m irrup 



(Z + s)!(Z-s)! 4tt 

I - s \ ( l + s \ , 
J ) 



1/2 / n \ 21 



sin 



j + s - m K ' \ 2 



2j+s— m 



(22) 



with the sum over j being understood in the range satisfying both < j < 
I — s and < j '• + s — m < I + s. In this limit, the eigenvalue A becomes 
A = A + O(auj) where A = 1(1 + 1) — s(s+l) is defined for later convenience. 

We may easily check that the radial equation (JT9|l reduces to the Teukol- 
sky equation j3HJ HOJ HIJ H2] when n = (hence /i = 2GM). The asymptotic 
solutions of eq. (|T9|) at the horizon and infinity are obtained in the same way 
as in four dimensions 



where 



r — > oo 


r - 


-> r h 


outgoing ingoing 


outgoing 


ingoing 






A~ s e~ ikr * 


k = to — 


ma 




r 2 + a 2 ' 





and r* is defined by r* — ► r for r — > oo and 



cZr* 
(ir 



r 2 + a 2 
A(r) ' 



(23) 



(24) 



(25) 



3.2 Hawking radiation and greybody factor 

Since we have shown that the massless brane field equations are separable 
into radial and angular parts, we may write down the power spectrum of the 
Hawking radiation [21] for each massless brane mode 



dE 



s,Lm 



dt du dip d cos i? 



1 s^i,mir h ,a]Lo) 

27Te ^-mn)/T _ 



,Sj m (ao>; i?, <p)\ 2 u}, (26) 



16 The so-called Condon-Shortley phase (— l)™ 1 is inserted to reduce into the standard 
definition of the spherical harmonics when s = 0: oYi m ($, ip) = Yj m ('(?, if). 



11 



where T and VL are the Hawking temperature and the angular velocity at the 
horizon, respectively given by 



T 



(n + l) + (n-l)q 2 
4ir{l + al)r h ' 



[l + al)r h 



(27) 



and s T^ m {rh,a;uj) is the greybody factor [23 EI] which is identical to the 
absorption probability of the incoming wave of the corresponding mode. (In 
this paper we only consider the modes which can be treated as massless 
compared with the Hawking temperature T since the emissions from massive 
modes are Boltzmann suppressed; Typically the standard model fields can 
be treated as massless at the LHC energy range.) Integrating eq. (j2B|) by the 
angular variables, we obtain 



dt du 



lm 



2ire( u - mfi )/ T - (-1) 



2.s 



U). 



(28) 



In the limit au <C 1 we can also write down the angular dependent power 
spectrum utilizing eq. (|2*T|) 



dE 



lm 



dt d cos d duj 



2 n e ("-mfi)/T _ 



2.s 



10 



2k 



d(p\ 8 Yi m (-&,(p)\ : 



(29) 



where the integral in the square brackets can be done with eq. (|22|): we 
summarize the results for the leading modes in the following table. 



s I m 


f^d<p\ s Y lm (^<p)\ 2 



1 1 
1 
1 -1 


1/2 
(3/4) sin 2 ■& 
(3/2) cos 2 
(3/4) sin 2 $ 


1/2 1/2 1/2 
1/2 1/2 -1/2 


sin 2 (^/2) 
cos 2 (tf/2) 


1 1 1 
1 1 
1 1 -1 


(3/8)(l -cos^) 2 

(3/4) sin 2 ■& 
(3/8) (1 + cost?) 2 



Approximately, the time dependence of M and J can be determined by 



d_ 

dt 



M 
J 



s,l,in 



du 



e (oj-mQ)/T _ ( — 1) 



0J 

m 



(30) 



where g s is the number of 'massless' degrees of freedom at temperature T, 
namely the number of degrees of freedom whose masses are smaller than T, 
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with spin s. (Typically g — 4, gi — 90 and g\ = 24 when T > m t ,m H 
and g — 0, gi — 78 and g\ = 18 when < T < mw for the standard 

2 

model fields.) Therefore, once we obtain the greybody factors, we completely 
determine the Hawking radiation and the subsequent evolution of the black 
hole up to the Planck phase, at which the semi-classical description by the 
Hawking radiation breaks down and a few quanta radiated is not predictable. 

In the high frequency limit, the absorption cross section for each mode a = 
(7r/co> 2 )r is supposed to reach the geometrical optics limit (see e.g. refs. [T31 

d) 

/ n + 3 x2/(n + l) n + 3 

°,o. = - [—) —r H . (31) 

In all the phenomenological literature this limit have been applied when 
one calculate the Hawking radiation. (In the refs. phenomenological 
weighting factors 2/3 and 1/4 are multiplied to eq. lJBTj) for s = 1/2 and s = 1 
fields, respectively, based on the result in four dimensions [HI].) 

3.3 Greybody factors for Randall- Sundrum black hole 

To obtain the greybody factors from eqs. ()18)) and ()19j) in general dimensions, 
we need the numerical calculation, which is beyond the scope of this paper 
and will be shown in ref . [S^ . In this paper we present analytic expression of 
the greybody factors of brane fields for n = 1 Randall- Sundrum black hole 
within the low frequency expansion. 17 Here we outline our procedure: First 
we obtain the "near horizon" and "far field" solutions in the corresponding 
limits; Then we match these two solutions at the "overwrapping region" in 
which both limits are consistently satisfied; Finally we impose the "purely 
ingoing" boundary condition at the near horizon side and then read the 
coefficients of "outgoing" and "ingoing" modes at the far field side; The ratio 
of these two coefficients can be translated into the absorption probability of 
the mode, which is nothing but the greybody factor itself. 
First for convenience, we define dimensionless quantities 

r — r h ~ uj — mVt 2 

£= , u = r h u, Q= = (1 + aju - ma*. (32) 

r h ZtyI 

(Note that in the Schwarzschild limit a* — > 0, Q becomes Q — > u5.) Then the 
radial equation (fTTJ|) becomes 

+ 2) 2 it!^ + 2(s + + !)(£ + 2)R^ + VR = 0, (33) 



17 



See ref. for the study of bulk scalar emission of five dimensional black hole. 
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where 



+ 2) + Q\ + 2isu£(£ + 1)(£ + 2) - 2zsQ(£ + 1) 
-[A + O(a*uj)m + 2). (34) 



V 



In the near horizon limit « 1, the potential ()34j) becomes 

V = Q 2 - 2is(£ + l)Q - A ^ + 2) + O(u>0, (35) 

and the solution of eq. with the potential (J35j) is obtained with the 
hypergeometric function 



R 



NH 



Ci 



T i4 



2 Fi(-Z - s, Z - s + 1, 1 - s - iQ; -|) 



2 Fi H + Z + 1 + iQ, 1 + s + iQ; 



To impose the ingoing boundary condition at the horizon (|23)1 . i.e. 

c?r* Q 



s „—ikr* 



k- 



(36) 



(37) 



we put C2 = and normalize Ci = 1 without loss of generality and then we 
obtain 



R 



NH 



r 



. iQ 



2j PxH - s, Z - s + 1, 1 - 8 - iQ; -|). 



(3* 



In the far field limit £ ^> 1 + \Q\, the eq. (J3Hj) becomes 
2(s + l} 



= R 



+ 



fl, (39) 
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and the solution is obtained by the Kummer's confluent hypergeometric func- 
tion 

R FF = B ie -^(^j 1 F 1 (/-s + l,2/ + 2;2^0 

+ B 2 e-& (|) xFi(-l - s, -21; 2iu>£), (40) 

where singularity from 21 being integer is regularized by the higher order 
terms in cD. 

Matching the NH and FF solutions (|38|) and (J4*U|) in the overlapping region 
1 + \Q\ < £ < 1/^, w e obtain 

B = T(2l + l)T(l-s-iQ) B = T(-2l - l)r(l - s - jg) 

1 r(/-s + i)r(/ + i-zQ)' 2 r(-z-s)r(-z-«Q) 

Then we extend the obtained FF solution toward the region £ 3> 1/65 

i?oo = r in e-^^ +y oute ^^ , (42) 

where 

r(2/ + i)r(2/ + 2) r(i- 3 -»g) 

r( / _ s + l)T(l + s + l)T(l + l-zQ) [ ] 

r(- 2 /)r(-2/-i) T(i- 3 -iQ) , l+s 

+ T(-l - s)T(-l + s) T(-l-iQ) [ ^ ' 

_ r(2Z + i)r(2/ + 2) r(i - s - ig) . ^ 
out ~ [r(/-s + i)F r(/ + i-.g) (4 ^ 

+ ^T^ 2 ^^ (43) 

Let us define R_ s as the solution of the equation obtained by the flip of 
the sign of s, i.e., s — > — s from eq. (|19j). When A rr = 2 as in n = 1 (or as 
in the limit r ^> r# in n > 2), we may obtain the conserved current in the 
same way as in the four-dimensional case 

J = A (R- s d r R* s - R*d r R- s ) + sA ir R„ s R* s , (44) 

which satisfies d r J = 0. In the limit r 3> r H , 

R s ~ Fine—V-^F^e^r- 25 - 1 , 

~ Z^e-^r- 1 + ^e^r 2 ^ 1 , (45) 
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where Z in = Y in \ s ^_ s and Z out = Y" ut| s ^- S , and J becomes 



J ~ 2iLu(Z m Y*-Z out Y; ut) 



(46) 



Therefore, we may calculate the greybody factor T (=the absorption proba- 
bility) in the same way as the Page's trick [3T] 



Y nu tZ. 



out " out 



Y m Z\, 



l-C 



1 + C 



(47) 



where 



C 



{4iu) 2l+1 ni + s )\(i-s)V 

4 I (20K2Z + 1)! , 



2i+l 



(48) 



with (a) n = Iln'=i( Q; + n> — 1) being the Pochhammer's symbol. 

For concreteness, we write down the explicit expansion of eq. (|4~7jl up to 
0(u) 6 ) terms 



oro,o 

0^1, m 



Au 2 - 8£ 4 + 0(ti e 



4Q£ S 

9^ 
16Q£ 5 



(l + Q 2 )+0(u5 6 ) ; 



2025 V 4 4 / V 



10 \ 



2 2 ' 



2 2 ' 



(l + 4Q 2 ) " y (l + 4Q 2 ) 2 + O(0 e 



36 I 9 9 ) V ' 



lP2,m 



i + g 2 ) +o(d; e 



i6Qur 

225 1 4 4 ) v ; 



(49) 



Note that subleading terms in £3 are already neglected when we obtain 
eq. (|47jl and the terms from these contributions are not written nor included 
in eqs. (|4*7j) and p5jl . We also note that the so-called s-wave dominance is 
maximally violated for spinor and vector fields since there are no / = modes 
for them. 
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3.4 Radiations from Randall- Sundrum black hole 



The greybody factor (|47jl is obtained from low-frequency expansions. In 
four dimensions, it is known that the greybody factors in the low-frequency 
expansion provide smaller value for the right hand side of eq. IJHOj) than the 
one from full numerical calculation [32J . Therefore in this paper we do not try 
to show the time evolution of the black hole nor the time- integrated result. 

In Figs. |2H3 we show the power spectrum (|2*8|) for spin 0, 1/2 and 1 
fields. The black lines are our results for a* = 0, 0.5, 1.0 and 1.5 utilizing the 
expression (|4Tjl with up to I < 7 modes, respectively from below to above at 
the left of the peak (and from above to below at the right of the peak). Note 
that our approximation is valid for the region satisfying both of a*oj < 1 and 
u) < 1, typically at the left of the peak. Two gray lines below and above 
are the corresponding power spectra in the geometrical optics limit (jHlj) with 
and without phenomenological weighting factor, respectively (2/3 for spinors 
or 1 /4 for vectors) [3 E] . 

In Figs. ISHTU1 we present the angular dependent power spectrum (J29|) for 
spin 0, 1/2 and 1 fields when a* = (a*) max = 1-5. The modes are taken up 
to I < 1. We observe that there is large angular dependence for spinors and 
vectors. Note that = corresponds the direction of the angular momentum 
of the black hole which is perpendicular to the beam axis. The angular 
dependence shown in Figs. IH^ 101 vanishes when we take the limit a* — > 0. 

4 Summary 

We have studied theoretical aspects of the rotating black hole production 
and evaporation. 

For production, we present an estimation of the geometrical cross section 
up to unknown mass and angular momentum loss in the balding phase. Our 
result of the maximum impact parameter 6 max is in good agreement with the 
numerical result by Yoshino and Nambu when the number of extra dimen- 
sions is n > 1 (i.e. within 6.5% when n = 1 and 1.5% when n > 2), though 
ours predicts same as the naive value in the Schwarzschild approximation 
frmax = Ts{M) when n = which is 24% larger than the numerical result. 
(Here we note that our refinement from the Schwarzschild approximation 
results in the enlargement of the production cross section, contrary to the 
previous claim in the literature.) Relying on this agreement, we obtain the 
(differential) cross section for a given mass and (an interval of) an angular 
momentum, which increases linearly with the angular momentum up to the 
cut-off value J max = 6 max M/2. This result shows that black holes tend to 
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be produced with large angular momenta. We also studied the possibility of 
the black ring formation and find that it would possibly form when there are 
many extra dimensions. 

For evaporation, we first calculate the brane field equations for general 
spin and for an arbitrary number of extra dimensions. We show that the 
equations are separable into radial and angular parts as the four-dimensional 
Teukolsky equations. From these equations, we obtain the greybody factors 
for brane fields with general spin for the five-dimensions (n = 1) Kerr black 
hole within the low-frequency expansion. We present the resultant power 
spectrum which is substantially different from the one with geometrical ap- 
proximation utilized in the literature. 

We address several phenomenological implications of our results. The 
production cross section of the black holes becomes larger than the one cal- 
culated from the Schwarzschild radius. The more precise determination of the 
radiation power is now available. We have shown that the black holes are pro- 
duced with large angular momenta and that the resultant radiations will have 
strong angular dependence for s = 1/2 and s = 1 modes which points per- 
pendicular to the beam axis while very small angular dependence is expected 
for scalar mode. When we average over opposite helicity states, the up-down 
asymmetry with respect to the angular momentum axis shown in Fig. El dis- 
appears (though there still remains angular dependence itself) jSH Ell ! We 
expect similar result for vector fields (which correspond to Fig. IIP)) . More 
quantitative estimation will need the greybody factors for arbitrary frequency 
calculated numerically. 
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A Separability of brane fields 



The various field equations in the four- dimensional Kerr background are 
known to be separable. This results from the special feature of the four- 
dimensional Kerr metric, that is, the vacuum metric which has a pair of 
degenerate principal null directions (Petrov type D). The four-dimensional 
metric considered in this paper is the induced metric of the totally geodesic 
probe brane in the higher-dimensional Kerr field. Though this brane metric 
turns out to be of Petrov type D, it is not the vacuum metric itself. Nev- 
ertheless, it happens that the massless fields on the brane are separable, as 
shown bellow. 

The induced metric on the three-brane in the (4 + n)-dimensional Kerr 
metric (with a single nonzero angular momentum) is given in terms of the 
Boyer-Lindquist coordinate system by 

g = (l - Ir 1 -) dt 2 + ^r 1 -" sin 2 #dtd<p - sin 2 (V 2 + a 2 + ^ |^ V 
A 

where 

£ = r 2 + a 2 cos 2 A = r 2 + a 2 - fir 1 '" 1 , (51) 

and the parameters \i and a are equivalent to the total mass M and the 
angular momentum J 

_ (n + 2)A n+2 n _ A n+2 fia 

M ~ IMG ' J ~ ~^G~- {b2) 

where A n+2 = 2^ n+3 ^ 2 /T((n + 3)/2) is the area of a unit (n + 2)-sphere. 

The direct calculation shows that the massless scalar field equation sep- 
arates on this background geometry. If we set tp = R(r)S('d)e~ lult+ ' imip , then 
V 2 y? — becomes 



1 d 




^sin + (a 2 u 2 cos 2 - m 2 esc 2 ■& + A) S = 0, (53) 



K 2 A 2r + (n- l)ar n ]K 

— + Aiur - i- — 

A A 



+ 2mauj — a uj — A 



R = 0, (54) 



where K = (r 2 + a 2 )u — am. We note that the Hamilton- Jacobi and massive 
scalar field equations are also separable though we do not shown them here; 
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a test particle on the brane has an additional conserved quantity (Carter 
constant) besides the energy and the angular momentum. 

To show the separability of higher spinor field equations, we work on the 
Newman- Penrose formalism [HI]. 18 We set the null tetrad as follows: 

n„ = S^-a sin 2 <&5$ - —8^ 



n' 



< - (r 2 + a 2 )5- - Xr^ ; 



r — ia cos d 



2V2( r + iacos^) L f v / Mj 2l/2 

= %. (55) 

where the bar denotes the complex conjugation. These are subject to the 
normalization: n^n'^ 1 = — m M m' M = 1, n M n M = n'n ,,x = m^m^ = n^m^ = 
n'^rn^ = 0. Alternative description is given by the two-component spinor 
o A , i A via the identification 

« n'", m", m'") <-> {p A o A \i A l A \ o A I A ' , i A o A ' ) , (56) 

with the symplectic structure e^ 45 = o^f 8 — i A o B , e 01 = eoi = 1. Each 
component of the spinor covariant derivative V^a' is denoted by 

(V o, V tf , V or , Vis) = (A D', 5, 5'), (5T) 

and the spin-coefficients are defined by 

D(o,l) = (eo — kl, — t'o — el), D'(o, l) = (— e'o — tl, —k'o + e't), 
S(o,i) = (f3o — <jl, —p'o — j3i), S'(o, i) — {—fl'o — pi, —cr'o + 0C). 

(58) 



Then, all the nonzero spin-coefficients are 
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ia sini) 1 pcoti? . iap 2 sin 

P= : a' P= — 7mri T 



2V2E ' r-mcostf' 2^ ' y/2 ' 

p = — , e=p-— A, r , p=T+p. (59) 



18 See e.g. ref. |58| for review of the Newman-Penrose formalism and spinor calculations. 
We follow the conventions of this reference. 

19 Though we have defined the spin-coefficients in spinor form, the tensor calculus would 
work better in actual computation. Sec eqs. (4.5.28) in ref. |58j for the equivalent tensorial 
definition for the spin-coefficient. 
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Here, let us consider the Weyl equation (s = 1/2) and the Maxwell equa- 
tion (s = 1) on this background brane metric. 

We define the component of the Weyl spinor ipA simply by -0o — ^ao a and 
■01 = ^a^ A - Then, each component of the Weyl equation \J AA ipA becomes 
explicitly 



Dip, - 5'iPo 



(e'-p')0o + (r-/?)^i. 



(60) 
(61) 



On the other hand, the Maxwell field is represented by the second-rank 
symmetric spinor 0ab, and its components are denoted by 0o 



I) abo a l b and 



1)abo a o b , 

^>abl A l E \ respectively Then, the source-free 



Maxwell equation V 4>ab — leads to 



' i -6'<t> = {2(3' - r')0 o + 2p0i - K(j) 2 

>2 - 8'<j>i = cr'0o - 2r'0i + (p - 2e)0 2 

= (P' - 2e')0O - + CX0 2 

£'0i-502 = -k'0o + 2p'0i + (2(3 - r)0 2 . 



(62) 
(63) 
(64) 
(65) 



The brane-induced metric turns out to be of Petrov type D, namely, the 
gravitational spinor \1/ abcd has only nonzero component \I/ 2 =\I/ abcdo a o b l c i D . 
Besides this condition, when n = a = K,' = a' = hold as in the present case, 
we have the identities for the differential operators 



[D - (p + l)e + e + qp - p](S - p(3 + qr) 

- [5 - (p + 1)0 + p'-f' + qr}(D - pe + qp) = 
\D' - (p + l)e' + e' + qp' - p'}(5' - p/3' + qr') 

- [5' - (p + 1)(3' + p-f + qr'} (D' - pe' + qp') 



0. 



for any pair of the numbers (p, q), where we have used the identities 

5D-D5 = (? -P' + (3)D + kD' -(p + e-e)5-a8', 
5'D' - D'S' = n'D + (f-P + f3')D' - a'5 - (p' + e' - e')8'. 



(66) 
(67) 



(68) 
(69) 



Applying (5 + (3' - f - t) to Eq. (jHH and (D + e - p - p) to Eq. ([FITT). 
subtracting one from the other, and using Eq. (JBTIj) for (p, g) = (—1, —1), we 
obtain the decoupled equation for ip 



(r 2 + a 2 ) 2 



A 



+ 



2r - 



(r 2 



— a 2 sin 2 § 
fa 2 ) A 



<9 2 



V 2 + a 2 



5 2 



1 



5 2 



A 



2A 



+ ia cos'd 



d_ 

dt 



+ 



dtdip 
i cos^ 



aA, 



2 A sin 2 



\ A sin 2 d ) dip 2 
d_ 

dip 
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1 d . cot 2 "d 1 n{n — l)nr n 1 

dr dr sin $ dti d$ 4 2 2 



>^o = 0. 
(70) 



If we set ip = R(r)S('d)e- iu)t+mv , then we obtain 



1 d ( . JS S 
I sin 1/ — 

sin ■& d$ \ d'd 



+ I a 2 J 1 cos 2 -& - 



m 



sin 2 $ 



auj cos d — 



mcosd 1 9 n 1 A 
. 2 - - - cot 2 + -+A]S = O, 

(71) 



A" 1 / 2 — f a 3/2 — \ 
dr \ dr J 



+ 



K 2 . i [2r + (n-l)/ir-"]ir N 
— + 2iur - - - 



n(n — l)/xr 



-n-1 



+ 2mato — a 2 cu 2 — A 



R = 0. 



(72) 



For the Maxwell field, applying (5 - (3 + f3' - f' - 2r) to Eq. flE2J and 
(Z) — e + e — 2p — p) to Eq. (jHSJ), subtracting one from the other, and using 
Eq. (jfifij) for (p, g) = (0, —2), we obtain 



(r 2 + a 2 ) 2 
A 



a 2 sin 2 d 



dt 2 



+ 



2a(r 2 + a 2 



A 



- 2a 



d 2 



+ 
+ 



/ir n [(n + l)r 2 + (n — l)a 2 ] 
A 



+ 2(r + ia cos^) 



<9t<9<£> 

Of 



+ 



A sin 2 



a[2r + (n — l)«r ™] 2icosi? 

sin 2 



A 



-A 2 



A*" ( > ^ + -t^ - 1 + n(n - l),r-^U = 0. 



(73) 



Set 0o = J R(r)5(^)e-^* +imv , then we have 



1 d ( . JS\ 
I sin 1/ — I 

sm<dd$\ d-d) 

( 9 9 9 „ ^ 2 „ 2m cost? , „ A „ 

+ a 2 u; 2 cos 2 1? - - 2au cos •& r^— cot 2 ■& + l + A\S = 0, 



sin i9 



sin 2 $ 



(74) 
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n(n — l)/ir n +2mauj — auj — A 



R = 0. 



(75) 



In summary, the spin-s massless field equation becomes 



1 d ( . n dS s 
I sin u — 

sin $ d$ \ d$ , 



+ 



[s — au cos $) 2 — (s cot "d + m esc $) 2 — s(s — 1) + A 



S = 0, 



(76) 



and 



+ 2mau — a 2 u 2 — A 



K 2 ( [2r + (n-l)fir- n ]K 

— — h s Aiujr — i n[n — l)ar 



R = 0. 
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Figure 1: Schematic picture for the condition of the black hole formation. 



0. 0.2 0.4 0.6 0.8 1. 

rh dE / dt dco 

3 • 10 " 3 
2 • 10 " 3 
1 • 10 " 3 


0. 0.2 0.4 0.6 0.8 1. 

rhCO 

Figure 2: Scalar (s = 0) power spectrum ThdE / dtdu vs. u = r^uj in lenear- 
linear plot. The gray line corresponds the geometrical optics limit. The black 
lines are our results for a* = 0, 0.5, 1.0 and 1.5 from above to below. Note 
that our approximation is valid for u < min(l, a^ 1 ). 
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Figure 3: Scalar (s = 0) power spectrum r^dE / dtdui vs. r h uj in log-log plot. 
See the caption of Fig. El for explantion. 
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Figure 4: Spinor (s = 1/2) power spectrum r^dE / 'dtdu vs. uj = r^LO in 
linear-linear plot. Two gray lines below and above correspond to the geo- 
metrical optics limit with and without the phenomenological weighting factor 
2/3, respectively. The black lines are our results for a* = 0, 0.5, 1.0 and 1.5, 
respectively from right to left at the peak location. Note that our approxi- 
mation is valid for cD < minfl, a! 1 ). 
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Figure 6: Vector (s = 1) power spectrum r h dE / dtdu vs. u = r h u in linear- 
linear plot. Two gray lines below and above correspond to the geometrical 
optics limit with and without the phenomenological weighting factor 1/4, 
respectively. The black lines are our results for a* = 0, 0.5, 1.0 and 1.5, 
respectively from below to above at the left of the peaks. Note that our 
approximation is valid for u < min(l, a^ 1 ). 
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Figure 7: Vector (s = 1) power spectrum r h dE / dtdu vs. r h uj in log-log plot. 
See the caption of Fig. El for explantion. 
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Figure 10: Vector (s = 1) power spectrum r^dE / dtdud cos $ vs. r^u and 
cost? for a* = 1.5. 
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